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1,  Introduction 

In  [1],  [2],  the  author  studied  the  regular  continued 
fraction  expansion 


(1.1)  [ ^2*  *•*  * 

for  the  binomial  quadratic  surd  ^C~  +  L ,  where  C  and  L  are  rational 

numbers,  C  positive.  Simultaneously,  the  approximations  x  to  n/C~+L 

n 

given  by  an  extension  of  Newton’s  formula, 


(1.2) 


x.  .  Xj  +  C  -  L2 
Xn  “xj  +Xj  -  2L 


where  xi  and  x.  are  certain  previous  approximations  to  the  value  of  Vc  +  L 
were  studied.  Let  A^/B^sx^^  denote  the  i-th  approximant  of  (1,1).  Then 
in  [  1],  [2],  a  complete  classification  was  derived  concerning  which  ones  of 
the  A./Bj  are  also  approximations  to  sfc~  +  L  given  by  formula  (1.2), 


Numbers  in  brackets  refer  to  the  bibliography  at  the  end  of  the  paper. 
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In  the  present  study,  the  use  of  automatic  digital  computers  is 
demonstrated  for  the  rapid  computation  of  the  regular  continued  fraction 
expansion  (1,1)  for  \Tc  +  L,  of  the  approximants  of  (1.1),  and  of  the 
approximations  by  formula  (1.  2)  .  The  use  of  an  electronic  computer  is 
especially  convenient  here  since  the  same  operations  are  repeated  many 
times  for  many  combinations  of  the  various  values  of  the  parameters. 

Simultaneously,  the  aim  of  this  paper  is  a  demonstration  of 
the  use  of  computers  in  problems  in  numerical  analysis,  and,  in  particular, 
in  their  application  to  continued  fractions  and  to  number  theory.  Very 
few  such  applications  have  thus  far  been  made.  Consequently,  they  are 
important  since  they  point  the  way  in  which  extensions  are  possible  in 
applications  of  continued  fraction  theory  in  functional  analysis,  and  in 
approximations  of  functions  by  rational  fractions  through  the  use  of 
continued  fractions. 

In  §2  are  given  the  mathematical  description  of  the  problem  to 
be  solved  and  the  formulas  and  methods  by  which  it  is  solved.  A  summary 
of  the  calculation  procedure  is  given  in  §3.  It  contains  an  Algol  program 
for  the  algorithms  involved  (Table  I)  and  a  sample  problem  computed  on 
the  PERM  machine  ( Technische  Hochschule,  Munchen)  ( Table  II) .  Of 
Interest  for  comparison  purposes,  Table  III  contains  a  flow  chart  and 
description  thereof  for  the  Fortran  program  of  the  previous  algorithms  ( Table  IV) . 
Table  V  shows  the  preceding  problem  worked  out  with  the  Fortran  program 
and  the  Control  Data  Corporation  1604  machine.  It  will  be  noted  that  the 
Fortran  program  is  ready  to  put  directly  into  the  machine.  On  the  other 


hand,  the  Algol  program  must  be  incorporated  into  a  main  program 
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2.  Mathematical  description 

I.  The  first  problem  is  the  derivation  of  formulas  for  the 
computation  of  the  numbers  b  .  in  the  regular  periodic  continued 
fraction  expansion  for  */c~+  L,  namely, 


(2.1) 


,  1  l 

0  bj+b^... 


ill 

+  b  +  b.  +  b  +  . . .  ’ 

p  i  2 


(formula  ( 1. 1) ) .  Lagrange  showed  that  this  continued  fraction  is 
always  periodic.  It  is  also  symmetric  if  L  =  0  ( cf. ,  for  example,  [  3] ) . 


SD+  P 

One  is  given  Vc~+  L.  This  can  be  written  ° 


Qo  ' 


where  (D  -  P  )  is  divisible  by  Q  .  Then  b  is  the  largest  integer 
o  o  o 


sTd+  P 

contained  in  - ^ — —  t  Thus,  one  is  given  the  values  of  D,  P  ,  Q  , 

O  -  o  o 


and  bQ.  The  continued  fraction  (2, 1)  is  generated  by  a  sequence  of 
linear  transformations  of  the  type 


(2.2) 


f  =  b  + 


v  v 


Now, 


fv  +  1  *  v  =  0  *  1 »  ••• 
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n/d"  +  P 


(2.3) 


v  h  ,  v+  1 

Qv  '  v  ^+Pv+l 


D  +  PP.+(P  +  P  =  bQ  */D  +  bQ  P  . ,  +  Q  Q 

V  V+l  V  V+l  V  V  V  V  V+l  v  v+l 


Since  *Td  is  irrational. 


(2.  4) 


D+P  P  =  b  Q  P  +  Q  Q  , 
v  v+l  v  v  v+l  v  v+l 


(2.5) 


P  +  P  =  b  Q  . 

V  V+l  V  V 


When  one  multiplies  equation  (2,5)  by  P^+1  and  subtracts  the  result 
from  ( 2. 4),  one  obtains 


D  -  Pv+1  *  QVQV+1  •  v=0-  1 . 


Then,  when  one  replaces  v  by  v  -  1  in  the  last  equation  and  subtracts 


the  result  from  (2.5),  one  obtains 


-.2  r.2 


QvQv+l  "  VlQv  -  Pv  '  Pv+1  =  (  Pv  '  Pv+1><  Pv  +  Pv+1)=  (Pv  -  "v+l)  bvQv  • 


(2.6) 


Q  . ,  -  Q  =  b  (P  -  P  ^, ), 

v+l  V-l  V  V  v+l 
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From  (2.5)  and  (2.6)  one  obtains  the  recurrence  relations 


(2.7) 


P  ,=  b  Q  -  P  , 

V+l  V  V  V 


'WW'WVl’'"0'1' 


Also 


\[D+  P 


(2.8) 


=  largest  integer  contained  in 


n/eT+  P 


Formulas  (  2.  7)  and  ( 2.  8)  are  those  required  for  the  computation  of 

the  P  ,  Q  ,  and  b  ,  v  =  1  ,  2 ,  ...  .  In  addition,  one  uses  the  recurrence 
v  v  v 

formulas  for  the  A  and  B  ,  the  numerator  and  denominator  of  the  v-th 
v  v 

approximant  of  (2.1): 


A  =  b  A  +  A  , B  =  b  B  ,  +  B 
v  v  v-1  v-2  v  v  v-1  v- 


2  * 


v  =  1  ,  2, 


(2.9) 


A  =  1  ,  A  -  b  ,  B  =  0  ,  B 
-1  oo  -l  o 


1 . 


It  is  remarked  that,  if  '/c"+  L  is  negative,  one  writes  -(-*/c~-  L) 
and  carries  out  the  above  computation  for  -*/C~ -  L.  The  minus  sign  is 
then  later  added  to  the  continued  fraction,  namely. 
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’  [  ^0  ’  bl  *  b2  '  ’**  ’  bp  ^ 


Furthermore,  it  is  understood  that  sTc  is  the  positive  root,  since, 


-  \Td+  P 


if  one  is  given 


,  one  writes  the  equivalent 


*fD  -  P 
_ o 

"Qo  ‘ 


Also  the  computation  involves  only  integers  if  one  first  multiplies 
numerator  and  denominator  of  the  binomial  quadratic  surd  in  question  by 


D  -  P 


such  a  constant  that 


=  Q  is  an  integer. 


In  order  to  find  the  period  p ,  one  notes  that  one  has  obtained 
a  complete  period  when  the  values  of  P^ ,  ,  and  b^  are  repeated. 

In  the  following  example  (Tables  Hand  V  ) ,  these  values  are  compared 

• 

with  P,  ,  Q.  ,  and  b.  since  there  is  only  one  initial  term  b  .  However, 
ill  o 

if  the  periodic  part  starts  with  b^  ,  j  =  2  ,  3 ,  ...  ,  one  uses  as  comparison 
numbers  P^ ,  ,  and  b^  . 

II.  Formula  (1.2)  is  simply  an  estimate  of  the  roots  of 
2  2 

x  -  2 Lx  -  (C  -  L  )  =  0.  It  is  an  extension  of  Newton's  formula.  In  [1] 
the  author  showed  that,  if  x  is  the  [(kp  -  s)j-th  approximant  of  (1.1) , 

Kp“  S 


x. .  .  .x  .  +  C  -  L 

x  _  (k-r)p-s  rp-1 _ 

kp-s  “  x  +  x  -  2L 


c, ,  .  +  x  , 

(k-r)p-s  rp-1 


(2.10) 
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k  =  2  ,  3,  ...  r  =  l,  2,  ...  k  -  r  >  1 ,  p  =  1 ,  2,  ..., 

where  x, ,  ,  and  x  .  are  the  [(  k  -  r)  p  -  s]-th  and  (rp-1)  -th 

(k-r)p-s  rp-1 

approximants,  respectively,  of  (1.1).  Sharper  formulas  for  xkp_g  , 

in  terms  of  i ,  were  obtained  in  [1]  and  [2]  when  (1.1)  is  symmetric 

of  even  (  p  =  2i)  and  odd  (  p  =  2i  +  1)  period.  However,  for  illustrative 

purposes  in  this  study,  formula  (2.10)  is  used. 

In  order  to  use  formula  (  2.10)  on  the  machine  to  compute 

\and  B  ,  one  must  write  the  formulas  for  these  two  quantities  in 
p-s  kp-s 

either  one  of  the  two  following  ways  (py  means  of  the  formulas  given  in  [1]  ) 


(2.11) 


\p-s  A(  k-r)  p-sArp-l  +  ^  bp  bo^  Arp-1  +  Arp  -2}  B(  k-r)  p-s 


Bkp-s  A(  k-r)  p-sBrp-l  +  ^p  bo^  Brp-1  +  Brp-2^B(  k-r)  p-s  ’ 


Q2A,  =QZA,,  v  A  . +  (D-P2)B,.  *  B  ., 

o  Kp-s  0  (k-r)  p-s  rp-1  o  (k-r)  p-s  rp-1 


(2.12) 


Q2B,  =Q2(A/t  .  B  .  +  A  ,B. .  .  )-2P  Q  B. ,  .  E 

vo  kp-s  0  (k-r)  p-s  rp-1  rp-1  (k-r)  p-s'  o  o  (k-r)  p-s 


rp-1  ' 


Similar  transformations  are  necessary  for  the  other  formulas  given  in 
[1]  and  [  2]  . 
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III.  As  an  estimate  of  the  error  involved,  the  formula  (cf.  (  1.3) 


of  [1]) 


(2.13) 


(\fc  +  L  ) 


2, 


9 


is  used  here. 


3.  Algol  and  Fortran  Programs  and  a  Sample  Problem. 


TABLE  I 

ALGOL  PROGRAM 

'PROCEDURE'  CONTINUED  FRACTION(  D,  Po,  Qo,  Co,  N)  RESULT:  ( A,  B,  C,  E,  PERIOD); 
'INTEGER'  D,  Po.Qo,  Co,  N,  PERIOD; 'INTEGER'  'ARRAY'  A,  B,  C;  'ARRAY'  E; 

'BEGIN'  'INTEGER'  NY;  'INTEGER'  'ARRAY'  P[o:N],  Q[  - 1 ;  N] ; 


P[o]  :  =  Po; 
Q[o]  :  =  Qo; 


A[  o]  :  =  C[  o]  :  =  Co; 


A[  -1]  ;  =  B[o]  ;  =  E[o]  ;  =1; 
B[  -1]  ;  =  PERIOD  :  =  o; 

Q[  - 1 J  ;  =  (D  -  PoxPo)/Qo; 
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'FOR'  NY  :  =  1  'STEP'  1  'UNTIL'  N  ’DO' 

'BEGIN' 

P[  NY]  :  =  C[NY  -  1]  xQ  [NY  -  1]  -  P[NY  -  1]  ; 

Q[NY]  :  =  C[NY  -  I]  x(  P[  NY  -  1]  -  P[  NY] )  +  Q[  NY  -  2]; 

C[  NY]  :  ENTIER((SORT(D)  +  P[  NY] ) /Q[  NY] )  ; 

A[  NY]  :  =  C[  NY]  xA[NY-l]  +  A[  NY  -  2]  ; 

B[  NY]  :  =  C[  NY]  xB  [NY  -  1]  +  B[  NY  -  2]  ; 

E[  NY]  :  =  l/(  B[  NY]  xB[  NY] )  ; 

'IF'  NY  'GREATER'  1  'AND'  (P[NY]  'EQUAL'  P[l]  'AND' 
Q[  NY]  'EQUAL'  Q[  1] )  'AND'  PERIOD  'EQUAL'  o  'THEN' 
PERIOD  :  =  NY  -  1 

'END' 

'END'  CONTINUED  FRACTION; 

•PROCEDURE'  EXTENSIONS,  Mf  A,  B,  C,  E,  PERIOD)  ; 

'INTEGER'  N,  M,  PERIOD;  'INTEGER'  'ARRAY' A, B,C;  'ARRAY'  E; 
'BEGIN'  'INTEGER'  K,S,  I,  NY; 

'IF1  PERIOD  'EQUAL'  0  'THEN' 

'BEGIN'  WRITE  (  "  NO  PERIOD  "  )  ;  'GOTO'  EXIT  EXTENSION  'END' 


'FOR'  NY  ;  =  N  +  1  'STEP'  1  'UNTIL'  M  'DO' 
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'BEGIN' 

K  :  =  ENTIER(  NY/PERIOD)  +  1  ; 

S  :  =  KxPERIOD-NY  ; 

1  :  =  (  K-I)  xPERIOD-S  ; 

A[  NY]  :  =  A[l]xA[  PERIOD  -  1]  +( (C[  PERIOD] -C[  o])  xA[  PERIOD  -  1]  + 
A[  PERIOD  -2])  XB[  1]  ; 

B[  NY]  :  =  A[  1 J  xB[  PERIOD  -  1 J  +  ((C[  PERIOD]  -  C[  o])  xB[PERIOD  -  1  j  + 
B[  PERIOD  -  2])  xB|  1]  ; 

E[  NY]  :  =  l/(  Bf  NY]  xB[  NY] ) 


'END'  ; 

EXIT  EXTENSION  :  'END'  EXTENSION  ; 

'PROCEDURE'  CORRECT  FORMULA ( D,  Po,  Qo,  M,  L,  A,  B,  C,  E,  PERIOD)  ; 

'INTEGER'  D,  Po,  Qo,  M,  L,  PERIOD;  'INTEGER'  'ARRAY'  A,  B,C;  'ARRAY'  E; 
'BEGIN'  'INTEGER'  K,  S,  I,  NY; 

'IF'  PERIOD  'EQUAL'  o  'THEN' 

'BEGIN1  WRITS  C  '  NO  PERIOD  1  '  ) ;  'GOTO'  EXIT  CORRECT  FORMULA 
'END'  ; 

'FOR'  NY  :  -  M  +  1  'STEP'  1  'UNTIL'  L  'DO' 


BEGIN 
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K  :  =  ENTIER{  NY/PERIOD)  +  1; 

S  :  =  KxPERIOD  -  NY; 

I  ;  =  (  K  -  1 )  xPERIOD  -  S; 

A[NY]  :  =  Z(l]xA[PERIOD-l]+((D-PoxPo)xB[l]xB[PERIOD-l  ])/(QoxQo)  ; 
B[  NY]  :  =  A[  1  ]  xB[  PERIOD  -  1]  +  A[  PERIOD  -  l]xB[  1]  -  (  2xPoxB[  1  ] 
xB[  PERIOD  -  1] )/ Qo; 

E[  NY]  :  =  l/{  B[  NY]  xB[  NY] ) 


’END*  ; 

EXIT  CORRECT  FORMULA  :  ’END*  CORRECT  FORMULA; 


Here  b  =  c  and  r  =  1 . 
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TABLE  II 

Algol  Program  and  PERM  Computer 


\T%  +  10 

4 


INPUTS: 

D:  96 

Po:  10 

Qo:  4  Co:  4 

SUBSCRIPT 

A 

B 

C 

l/( BxB) 

RESULTS  OF  PROCEDURE 

CONTINUED  FRACTION: 

0 

4 

1 

4 

+1.00000000 

1 

5 

1 

1 

+1.00000000 

2 

94 

19 

18 

+0.00277008 

3 

99 

20 

1 

+0.00250000 

4 

391 

79 

3 

+0.00016023 

5 

490 

99 

1 

+0.00010203 

6 

9211 

1861 

18 

+0.00000029 

7 

9701 

1960 

1 

+0.00000026 

8 

38314 

7741 

3 

+0.00000002 

9 

48015 

9701 

1 

+0.00000001 

10 

902584 

182359 

18 

+0.00000000 

RESULTS  OF  PROCEDURE 

EXTENSION: 

11 

950599 

192060 

+0.00000000 

12 

3754381 

758539 

+0.00000000 

RESULTS  OF  PROCEDURE 

CORRECT  FORMULA: 

13 

4704980 

950599 

+0.00000000 

14 

88444021 

17869321 

+0.00000000 

15 

93149001 

18819920 

+0.00000000 
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Description  of  the  Flow  Chart 

Box 

1  The  index  i  is  set  equal  to  zero. 

2  The  quantities  Q.  ,  D,  A^  ,  A^  B._1  ,  B.  are  read  as  input. 

is  computed  and  designated  as  b  . 

i 


3  The  quantity 


VF+  P. 

Q~ 


4  The  index  i  is  increased  by  1 . 

5  The  quantity  b^Q  i  .  i  ‘  Pi  .  i  is  computed  and  from  here  on  is 
designated  as  P^ . 

6  If  i  =  1  ,  box  7  is  carried  out;  if  i  *  1  ,  box  8  is  carried  out. 

2 

D-Pi 

7  The  quantity  -rr -  is  computed  and  from  here  on  is  designated  as  Q  . 

Wi-1  1 


8  The  quantity  b^  P^  -  P^)  +  2  is  computed  and  from  here  on  is 


designated  as  Q^. 


9  a)  The  quantity 


as  bj. 


4d+  p 


is  computed  and  from  here  on  is  designated 


b)  The  quantity  b  A  +  A  is  computed  and  from  here  on  is  designated 
1  *  • 


as  Ai# 


c)  The  quantity  b^^  +  B1  2  is  computed  and  from  here  on  is  designated 


as  Bi. 


#363 


-15- 


10  If  i  =  n,  the  instruction  in  box  11  is  carried  out;  if  i  *  n,  the 
computation  begins  again  with  box  4. 

11  The  results  are  printed. 

12  The  value  is  set  equal  to  2. 

13  If  P  =  P  ,  Q  =  Q  ,  and  b  =  b  ,  the  computation  proceeds 

1  i  1 

to  box  15;  otherwise,  the  computation  proceeds  to  box  14. 

14  The  quantity  is  increased  by  1 . 

15  a)  The  quantity  -  1  is  computed  and  from  here  on  is  designated  as  p. 

b)  The  quantity  p  -  1  is  computed  and  from  here  on  is  designated  as  s . 

c)  The  quantity  r  is  set  equal  to  1  . 

d)  The  quantity  K  is  set  equal  to  2 , 

e)  The  quantity  (K  -  r)p  -  s  is  computed  and  from  here  on  is  designated  as  i1  . 

f)  The  quantity  rp  -  1  is  computed  and  from  here  on  is  designated  as  i^ . 

g)  The  quantity  Kp  -  s  is  computed  and  from  here  on  is  designated  as  1, 

2  2  2 

h)  The  quantity  (Q  (A  A.  )  +  ( D  -  P  )  B.  B.)/Q  is  computed  and  from 

o  iji2  o  ij  i2  o 

here  on  is  designated  as 

j 

2  2 

i)  The  quantity  ( Q  ( A.  B.  +  A.  B.  )  -  2P  Q  B.  B.  )  /Q  is  computed 

°  *1  z  iz  h  0  0  il  iz  ° 

and  from  here  on  is  designated  as  B^,. 


16  The  results  are  printed  and  the  computation  ended, 
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TABLE  IV 

FORTRAN  PROGRAM 

THE  EXPRESSION  OF  A  QUADRATIC  SURD  AS  A  CONTINUED  FRACTION. 
PROGRAM  DONE  ON  CDC  1604 

DIMENSION  LP(  50)  ,  LQ(  50)  ,  LB(  50)  ,  LA(  50)  ,  LC(  50) 

LLL  =  5 

IA  =  14073000 

IA  a  IA*10000000 

1  READ  2  ,  LPZ  f  LQZ ,  LD ,  N 

2  FORMAT  (3110,  14) 

200  FORMAT  (  1H  16  ,  116,  4121) 

PRINT  3,  LPZ,  LQZ,  LD,  N 

3  FORMAT  (4H1PO  =110,  8H  QO=I10,7H  D  =  I10,7H  N  =  14) 
M  =  SQRTF  ( FLOATF  (  LD ) ) 

5  LCZ  =  ( LPZ  +  M)  /LQZ 
7  LP  ( 1)  =  LCZ*LQZ  -  LPZ 

LQ  (1)  =  (LD  -  LP(  1)  *LP(  1)) /LQZ 
IF  (LQ(1))10,  28,  10 

10  LC  ( 1)  =  (  M  +  LP(  1) )  /LQ(  1) 

11  LA  ( 1)  *  LC  ( 1)  *  LCZ  +  1 
LAZ  =  LCZ 

LBZ  =  1 

LB  ( 1)  =  LC  ( 1) 

LP(  2)  =  LC  ( 1)  *  LQ(  1)  -  LP  ( 1) 
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LQ(2)  a  LC(  1)  *{LP(  1)  -  LP(  2))  +  LQZ 
IF(LQ(  2))  14,  28,  14 

14  LC(  2)  =  (  M  +  LP(  2))/LQ{  2) 

15  LA(2)  =  LC(2)  *LA(  1)  +  LAZ 

LB(  2)  =  LC(  2)  *LB(  1)  +  LBZ 

I  =  0 

PRINT  100 

100  FORMAT  (116H  SUBSCRIPT  P  Q 

1  C  A  B  ) 

PRINT  200,  I,  LPZ,  LQZ,  LCZ,  LAZ,  LBZ 
IB  =  1 

PRINT  200,  IB,  LP(  1 ) ,  LQ{  1) ,  LC(  1)  ,  LA(1),  LB(  1) 

IC  a  2 

PRINT  200,  IC,  LP(  2) ,  LQ(  2) ,  LC(2),  LA(  2) ,  LB(  2) 

DO  22  J  =  3,  N 

LP(J)  =  LC(  J  -  1)*LQ(J  -  1)  -  LP(J  -  1) 

LQ(J)  3  LC(J  -  1)*(LP(J  -  1)  -  LP(J))  +  LQ(  J  -  2) 

IF(  LQ(  J))  19,  28,  19 

19  LC(J)  =  (M+  LP(J))/LQ(J) 

20  LA( J)  =  LC(J)*LA(J  -  1)  +  LA(J  -  2) 

IF{  LA(  J)  -  IA)  21,  21,  29 

21  LB(  J)  a  LC(J)  *LB(J  -  1)  +  LB(J  -  2) 

IF(  LB( I)  -  IA)  30,  30,  31 

30  IC  a  IC  +  1 


PRINT  200,  IC,  LP(  J) ,  LQ(  J) ,  LC(J),  LA(J),  LB(J) 
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22  CONTINUE 

28  PRINT  800 

800  FORMAT{  7H0Q(  J)  =  0) 

GO  TO  190 

29  PRINT  900 

900  FORMAT(  14H  OVERFLOW  ON  A) 

GO  TO  190 

31  PRINT  901 

901  FORMAT  ( 14H  OVERFLOW  ON  B) 

GO  TO  190 

190  PRINT  191 

191  FORMAT  (43H  MORE  A  AND  B  USING  NEWTON  EXTENDED  FORMULA) 
PRINT  192 

192  FORMAT  (31H  SUBSCRIPT  A  B) 

KK  =  2 

599  IF ( LP(  1 )  -LP(KK))  600,  601,  600 

600  KK  =  KK  +  1 
GO  TO  599 

601  IF(  LQ(  1)  -  LQ{  KK) )  602,  603,  602 

602  KK  =  KK  +  1 
GO  TO  599 

603  IF{  LC(  1)  -  LC{  KK) )  604,  605,  604 

604  KK  =  KK  +  1 


GO  TO  599 


605  MP  =  KK  -  1 

MS  =  MP  -1 
IF(  MS)  505,501,505 
505  MR  =  1 

MK  =  2 

DO  501  IS  =  1  ,  MS 

11  =  (  MK  -  MR)  *  MP  -  IS 

12  =  MR*  MP  -  1 
11  =  MK  *  MP  -  IS 

LA(II)  =  (LQZ**2)*LA(I1  )*LA(I2)+  ( LD  -  LPZ**2)*LB(  II  )*LB(  12 ) 

LB(II)  =LQZ**2*(LA(I1  )*LB(  I2)  +  LA(  12)  *LB(  1 1))  -  2*LPZ*LQZ*LB(  II  )*LB(I2 
1) 

LA(II)  =  LA(  II)/LQZ**2 
LB(  II)  =  LB{  II)/LQZ**2 
PRINT  333,  II,  LA(  II),  LB(II) 

333  FORMAT<I8,  2116) 

501  CONTINUE 

IF{  LLL) 9 1 1 »  700,  700 
700  LLL  =  LLL  -  1 

GO  TO  1 
911  STOP 

END 


END 


-20- 


9  363 


TABLE  V 


Fortran  Program  and  Control  Data  1604 


PO  = 

10 

QO  = 

4 

D  =  96  N  =  30 

SUBSCRIPT 

P 

Q 

C 

A 

B 

0 

10 

4 

4 

4 

1 

1 

6 

15 

1 

5 

1 

2 

9 

1 

18 

94 

19 

3 

9 

15 

1 

99 

20 

4 

6 

4 

3 

391 

79 

5 

6 

15 

1 

490 

99 

6 

9 

1 

18 

9211 

1861 

7 

9 

15 

1 

9701 

1960 

8 

6 

4 

3 

38314 

7741 

9 

6 

15 

1 

48015 

9701 

10 

9 

1 

18 

902584 

182359 

11 

9 

15 

1 

950599 

192060 

12 

6 

4 

3 

3754381 

758539 

13 

6 

15 

1 

4704980 

950599 

14 

9 

1 

18 

88444021 

17869321 

15 

9 

15 

1 

93149001 

18819920 

16 

6 

4 

3 

367891024 

74329081 

17 

6 

15 

1 

461040025 

93149001 

18 

9 

1 

18 

8666611474 

1751011099 

19 

9 

15 

1 

9127651499 

1844160100 

20 

6 

4 

3 

36049565971 

7283491399 

21 

6 

15 

1 

45177217470 

9127651499 

22 

9 

1 

18 

849239480431 

171581218381 

23 

9 

15 

1 

894416697901 

180708869880 

24 

6 

4 

3 

3532489574134 

713707828021 

25 

6 

15 

1 

4426906272035 

894416697901 

26 

9 

1 

18 

83216802470764 

16813208390239 

27 

9 

15 

1 

87643708742799 

17707625088140 

28 

6 

4 

3 

64672951988506 

69936083654659 

OVERFLOW  ON  A 

MORE  A  AND  B  USING  NEWTON  EXTENDED  FORMULA 

SUBSCRIPT  A  B 

9701  I960 

9211  1861 

490  99 


7 

6 

5 
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